Ionization of a Model Atom: Exact Results and Connection with Experiment 
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We prove that a model atom having one bound state will be fully ionized by a time periodic 
potential of arbitrary strength r and frequency uj. The survival probability is for small r given by 



for times of order T 



where n is the number of "photons" required for ionization, 



with enhanced stability at resonances. For late times the decay is like t . Results are for a Id 
system with a delta function potential of strength —g(l + rj(t)) but comparison with experiments on 
the microwave ionization of excited hydrogen atoms and with recent analytical work indicate that 
many features are universal. 

PACS: 32.80.Rm, 03.65.Db, 32.80.Wr. 
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Transitions between bound and free states of a system 
are of great importance in many areas of science Q and 
"much of the practical business of quantum mechanics 
is calculating exponential decay rates" There are, 
however, still many unresolved questions when one goes 
beyond perturbation theory 0. Unfortunately, ap- 
proaches going beyond perturbation theory such as Flo- 
quet theory, semi-classical analysis and numerical solu- 
tion of the time dependent Schrodinger equation are both 
complicated and also involve, when calculating transi- 
tions to the continuum, uncontrolled approximations jjj]- 
||. It is only recently that some general results going 
beyond perturbation theory have been rigorously estab- 
lished for models with spatial structure . We still don't 
know, however, many basic facts about the ionization 
process, e.g. the conditions for a time dependent exter- 
nal field to fully dissociate a molecule or ionize an atom, 
much less the ionization probability as a function of time 
and of the form of such a field || . Granted that the prob- 
lem is intrinsically complicated it would be very valuable 
to have some simple solvable models which contain the 
spatial structure of the bound state and the continuum 
and can thus serve as a guide to the essential features of 
the process. 

In this note we describe new exact results relating to 
ionization of a very simple model atom by an oscillat- 
ing field (potential) of arbitrary strength and frequency. 
While our results hold for arbitrary strength perturba- 
tions, the predictions are particularly explicit and sharp 
in the case where the strength of the oscillating field is 
small relative to the binding potential — a situation com- 
monly encountered in practice. Going beyond perturba- 
tion theory we rigorously prove the existence of a well 
defined exponential decay regime which is followed, for 
late times when the survival probability is already very 
low, by a power law decay. This is true no matter how 
small the frequency. The times required for ionization 
are however very dependent on the perturbing frequency. 
For a harmonic perturbation with frequency uj the log- 
arithm of the ionization time grows like r~ 2n , where r 
is the normalized strength of the perturbation and n is 



the number of "photons" required for ionization. This 
is consistent with conclusions drawn from perturbation 
theory and other methods (the approach in || being the 
closest to ours), but is, as far as we know, the first exact 
result in this direction. We also obtain, via controlled 
schemes, such as continued fractions and convergent se- 
ries expansions, results for strong perturbing potentials. 

Quite surprisingly our results reproduce many features 
of the experimental curves for the multiphoton ionization 
of excited hydrogen atoms by a microwave field || . These 
features include both the general dependence of the ion- 
ization probabilities on field strength as well as the in- 
crease in the life time of the bound state when —nhuj, n 
integer, is very close to the binding energy. Such "reso- 
nance stabilization" is a striking feature of the Rydberg 
level ionization curves [3]. These successes and compar- 
isons with analytical results 0- J^] suggest that the sim- 
ple model we shall now describe contains many of the 
essential ingredients of the ionization process in real sys- 
tems. 

The model we consider is the much studied one- 
dimensional system with Hamiltonian || , Q , Q , 

h 2 d 2 

H = - — — -gd{y), g > 0, -oo < y < oo. (1) 
Am dy z 



Hq has a single bound state Ub(y) — ^/poe 
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Po = 



j^g with energy —tkoo = —h pg/2m and a continuous 
uniform spectrum on the positive real line, with general- 
ized eigenfunctions 



u(k,y) 
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-oo < k < oo 



Pa + i\k\ 

and energies h 2 k 2 /2m. 

Beginning at some initial time, say t = 0, we apply a 
perturbing potential —gr)(t)5(y), i.e. we change the pa- 
rameter g in Hq to g(l + r)(t)) and solve the time depen- 
dent Schrodinger equation for tp(y,t), 

v J {y,t) = 6{t)u b {yy^ t 

r°° 2 
+ / e(k,t)u(k,y)e-^ t dk (t>0) (2) 
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with initial values 0(0) = 1, 9(fc,0) = 0. This gives 
the survival probability |6>(i)| 2 , as well as the fraction of 
ejected electrons \Q(k,t)\ 2 dk with (quasi-) momentum in 
the interval dk. 

In a previous work H we found that this problem can 
be reduced to the solution of a single integral equation. 
Using units in which po,ujq, h, 2m and § equal 1 we get 



6{t) = 1 + 2i J,,' Y(s)ds 
Q(k, t) = 2|fe|/[v^F(l - i\k\)] Jo Y( s )e^ 1+k2 

where Y(t) satisfies the integral equation 



"ds 



Y(t) = V{t) { 1 + J [2i + M(t - t')]Y{t')dt' 



(3) 
(4) 



(5) 



To understand the ionization processes as a function of 
t, uj, and r requires a detailed study of the singularities 
of y{p) in the whole complex p-plane. This yields the 
following results: For small r, yip) has square root branch 
points at p = {— i(ncj + 1) + (3(r 2 ) : n E Z}, is analytic 
in the right half plane and also in an open neighborhood 
M of the imaginary axis with cuts through the branch 
points. As \q\ — > oo in A^ we have \y(q)\ = Oiru>\q\~ 2 ). 

If \u>— ^1 > const. r 2 , n a positive integer, then for small 
r the function y is meromorphic in the strips —muj — 1 — 
0(r 2 ) > 3(p) > -tow - w - 1 + 0(r 2 ), m e Z and has 
a unique pole in each of these strips, at a point p with 
> K(p) = 0(r 2 ") for small r. JS iften /oZ/ows f/iaf 0(i) 
can be decomposed as pi 



with 



Bit) = e-rt^WFuft) + 
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M(s) = - 

7T ./0 



du=-\l- 



,3/2 



An important result of the present work is that when 
rjit) is a trigonometric polynomial with real coefficients 



r)(t) = Aj smijujt) + Bj cos(jwt) 



(6) 



tends to zero as t 



the survival probability |$(£)| 2 
/or all lo > 0. 

This result follows from ||) and (H) once we establish 
that 2|Y(i)| = \0'(t)\ — > in an integrable way, and this 
represents the difficult part of the proof. Since the main 
features of the behavior of yip) are already present in 
the simplest case ij — r sin(wt) we now specialize to this 
case. The asymptotic characterization of Y is obtained 
from its Laplace transform yip) = J e~ pt Yit)dt, which 
satisfies the functional equation (cf. ([5])) 



yip) 



yip- 



yip- 



^/l — ip + uj — 1 



a/1 — ip — lu — 1 
no 



(7) 
w + p 

with the boundary condition yip) — > as 3(p) — » ±oo 
(the relevant branch of the square root is (1 — ip— w) 1 / 2 = 
— 1 + ip) 1 ! 2 for u> > 1). We show that the solution 
of (p]) with the given boundary conditions is unique and 
analytic for 5R(p) > 0, and its only singularities on the 
imaginary axis are square-root branch points (see below). 
This in turn implies that |V(£)| does indeed decay in an 
integrable way. The proof depends in a crucial way on the 
behavior of the solutions of the homogeneous equation 
associated to (0): y(p) has poles on a vertical line if the 
homogeneous equation has a solution that is uniformly 
bounded along that line. The absence of such solutions 
in the closed right half plane is shown by exploiting the 
symmetry with respect to complex conjugation of the 
underlying physical problem and carries through directly 
to the more general periodic potential (6). 



where F^ is periodic of period 2ttoj~ 1 
coefficients decay faster than r n n~ n / 2 , 
const.r m t~^l 2 for large t uniformly in m 



and its Fourier 
and \h m (t)\ < 
Furthermore, 
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for large t. 
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FIG. 1. Plot of log 10 |#(£)| 2 vs. time in units of ui^ 1 for 
several values of uj and r. Inset shows detail of power-law tail 
for uj = 1.3, r = 0.3. 

Consequently, for times of order l/3ff ('y) the survival 
probability decays as exp(— Tt), T — 25R(7), after which 
its long time behavior is |#(i)| 2 = 0(i -3 ). This is il- 
lustrated in Figure 1 where it is seen that for small r 
exponential decay holds up to times at which the sur- 
vival probability is extremely small, after which |#(i)| 2 
decays polynomially with many oscillations. Note that 
even for r as large as 0.3 the decay is essentially purely 
exponential for all practical purposes. Thus, for lo > 1 
Fermi's golden rule works magnificently B. 
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Using a continued fraction representation of the solu- 
tions of the homogeneous equation associated to (0) we 
obtain as r — > 0, 



r = <^ 



LO 

V2w - 1 



(1 - VT^Zj) 2 8w' 



2- 2n+2 y/nuj - 1 r 2 



ifw> l + 0(r 2 ) 
if we (i,l)+ 

if w e ( i ,^-r) + 



(9) 

where lo € (a, 6) + means a + 0(r 2 ) < < 6 — (^(r 2 ). The 



result for lo > 1 agrees with perturbation theory |lj] since 
the the transition matrix element is 



r 2 ). 1 
§sh 



< u b {y)\5{y)u{k,y) > 



1 



2vr 1 



(10) 



In Figure 2 we plot the behavior of T _1 which is just 
the time needed for |#(i)| 2 to decay significantly, as a 
function of lo. 




FIG. 2. log 10 F as a function of ui/uo at r = 0.01. At 
r = 0.2 (inset) shift of the resonance is visible. 

The curve is made up of smooth (roughly self-similar) 
pieces for lo in the intervals (n _1 , (n — l) -1 ) correspond- 
ing to ionization by n photons. Note that at resonances, 
when w _1 is an integer (i.e. multiple of lOq 1 here set equal 
to unity), the coefficient of r 2 ™, the leading term in T, 
goes to zero. At such values of lo one has to go to higher 
order in r, corresponding to letting lo approach the res- 
onance from below. This yields an enhanced stability of 
the bound state against ionization by perturbations with 
such frequencies. The origin of this behavior is, in our 
model, the vanishing of the matrix element in (|l0|) at 
k = 0. This behavior should hold quite generally since 
the quasi-free wavefunction u(k,y) may be expected to 
vanish pointwisc as k — > 0. For d > 1 there is an ad- 
ditional factor k d ~ 2 coming from the energy density of 
states near fc = 0. As r increases these resonances shift 



in the direction of increased frequency. For small r and 
lo = 1 the shift in the position of the resonance, some- 
times called the dynamic Stark effect [1], is about 

In Figure 3 we plot the strength of the perturbation 
r, required to make |#(i)| 2 = \ for a fixed number of os- 
cillations of the perturbing field (time measured in units 
of lo^ 1 ) as a function of lo. Also included in this figure 
are experimental results for the ionization of a hydrogen 
atom by a microwave field. In these still ongoing beau- 
tiful series of experiments, carried out by several groups 
and reviewed in ||, the atom is initially in an excited 
state with principal quantum number no ranging from 
32 to 90. The experimental results in Fig. 3 are taken 
from Table 1 in ||, see also Figures 13 and 18 there. The 
"natural frequency" loq is there taken to be that of a tran- 



sition from no to no + 1, loq 



The strength of the 



microwave field F is then normalized to the strength of 
the nuclear field in the initial state, which scales like Hq -4 . 
The plot there is thus of n^F vs. UqLo. To compare the 
results of our model with the experimental ones we had 
to relate r to UqF. Given the difference between the hy- 
drogen atom Hamiltonian with potential Vo(R) = — 1/R 
perturbed by a polarized electric field V\ = xF sm(Lot), 
and our model with V\ — rVo, this is clearly not some- 
thing that can be done in any unique way. We therefore 
simply tried to find a correspondence between n^F and 
r which would give the best visual fit. Somewhat to our 
surprise these fits for different values of lo/loq all turned 
out to have values of r close to Sn^F. A correspondence 
of the same order of magnitude is obtained by comparing 
the perturbation- induced shifts of bound state energies in 
our model and in Hydrogen. 
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FIG. 3. Comparison of the theoretical and experimental 
threshold amplitudes for 50% ionization vs. lo/ljo. 

The shift in the position of the resonances from the 
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integer fractional values seen in Fig. 2, due to the finite 
value of r, was approximated in Fig. 3 using the average 
value of r over the range, r w 0.195. 

In Figure 4 we plot |6*(£)| 2 vs. r for a fixed t and two 
different values of u). These frequencies are chosen to cor- 
respond to the values of lu/loq in the experimental curves. 
Figure 1 in |11| and Figure lb in Q. The agreement is 
very good for w/loq w .1116 and reasonable for the larger 
ratio. Our model essentially predicts that when the fields 
are not too strong, the experimental survival curves for a 
fixed 71qU> (away from the resonances) should behave es- 



sentially like exp ^— C[n^F] n 'o^ tcu j with C depending 

on UqUJ but, to first approximation, independent of n^F. 

The degree of agreement between the behavior of what 
might be considered as the absolutely simplest quantum 
mechanical model of a bound state coupled to the con- 
tinuum and experiments on hydrogen atoms is truly sur- 
prising. The experimental results and in particular the 
resonances have often been interpreted in terms of clas- 
sical phase space orbits in which resonance stabilization 
is due to KAM-like stability islands ||. Such classical 
analogs are absent in our model as in fact are "photons" . 
On the other hand, the special nature of the edge of the 
continuum seems to be quite general, cf. H. 



/ eo/(o =.434 
' r/(njF)=2.66 



(o/(i) =.1116 
r/(nSF)=3.85 
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FIG. 4. Ionization fraction at fixed t (corresponding to 300 
oscillations) as a function of amplitude. 

We note that for lu > luq, in the limit of small am- 
plitudes r, a predominantly exponential decay of the 
survival probability followed by a power-law decay was 
proven in j?J for three dimensional models with quite gen- 
eral local binding potentials having one bound state, per- 
turbed by a local potential of the form r cos(u)t)Vi (y). It 
seems likely that our results for general uj and r, including 
general periodic (perhaps also quasi-periodic) perturba- 



tions would extend to a similarly general setting. We are 
currently investigating various extensions of our model 
to understand the effect of the restriction to one bound 
state. This will hopefully lead to a more detailed under- 
standing, and some control over the ionization process. 

Because T relates to the position of the poles of the 
solution of (Q), a convenient way to determine T (mathe- 
matical rigor aside), if r is not too large, is the following, 
see also ||. One iterates n times the functional equa- 
tion (Q), n appropriately large, to express y(p) only in 
terms of y(p ± miuj) with \m\ > n. After neglecting the 
small contributions of the y(p ± mito), the poles of y(p) 
can be obtained by a rapidly converging power series in 
r, whose coefficients are relatively easy to find using a 
symbolic language program, although a careful monitor- 
ing of the square-root branches is required. A complete 
study of the poles and branch-points of y leads to (||) 
which is effectively the Borel summation of the formal 
(exponential) asymptotic expansion of Y for t — > oo. 
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